Uncovering wind turbine properties through two-dimensional stochastic modeling of wind dynamics 
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Using a method for stochastic data analysis, borrowed from statistical physics, we analyze synthetic data from 
a Markov chain model that reproduces measurements of wind speed and power production in a wind park in 
Portugal. From the theoretical point of view we argue that our methods can be used to extract unknown func- 
tional relations between two variables. We first show that indeed our analysis retrieves the power performance 
curve, which yields the relationship between wind speed and power production and discuss how such procedure 
can be extended for extracting functional relationships between pairs of physical variables in general. Second, 
we show how specific features, such as the turbine rated wind speed or the descriptive wind speed statistics, 
can be related with the equations describing the evolution of power production and wind speed at single wind 
turbines. 
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I. INTRODUCTION 

The use of efficient and clean renewable energy sources is 
one of the major conditions required to achieve the important 
aim of sustainable development in modern societies 1 1 1. Wind 
energy is one of such sources and wind turbines are being sub- 
ject to intensive studies for improving their efficiency 121 [3]. 
Although the basic laws of atmospheric wind motion have 
been known for a long time, important problems such as tur- 
bulence, layering, and the statistics of extreme events remain 
poorly understood. A better understanding of these phenom- 
ena can help to construct energy convertion schemes that are 
both more efficient and robust. Here, robustness must be con- 
sidered under two aspects: first, the occurrence of sudden 
changes in wind speed and direction can interrupt the pro- 
cess of energy conversion, meaning unreliability and a sud- 
den slump in the electrical energy generated, which is seen 
as one of the major obstacles for the replacement of fossil 
and nuclear plants by wind energy sources. Second, these 
sudden changes introduce massive mechanical stresses which 
can lead to excessive wear or, ultimately, to the destruction of 
wind generators. 

Wind flow is in general turbulent and non-homogeneous 21 
with a non-negligible stochastic contribution. Therefore, in 
order to be able to construct more accurate models for its 
physical properties, one needs either accurate measurements 
of the wind speed on the length scales of wind turbines, or 
suitable models that can statistically reproduce these mea- 
sured data. Since the wind turbines are driven by turbulent 
wind fields, the stochasticity of the wind fields transfers to 
stochastic dynamics of the wind turbine as a whole, of the 
loads on its structures and, last but not least, of the power 
output. Recently, a Markov chain model was used to repro- 
duce wind measured dataQ, based on the transition matrix 
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FIG. 1: Illustration of Iberia Peninsula indicating the position of the 
Portuguese wind park. In the inset one sees (bullets) the geographic 
location of each of the 57 wind turbines. In blue (AOl) one sees the 
wind turbine analyzed here. 



and time propagators for the wind speed and direction to- 
gether with the power production. Differently from previous 
first-order approaches|6J, information from two and three step 
transition probabilities are considered. 

In this paper, we aim at understanding the stochastic as- 
pects of power production coupled to the wind velocity field. 
To that end, we use a methodology recently introduced by 
some of us|7|, for uncovering optimal stochastic variables 
weakly|8| and strongly coupled[7J, and adapt it with two 
purposes. First, to properly derive the functional relation of 
pairs of variables whose values are extracted from the Markov 
chain model for wind turbines. Second, to uncover specific 
features of the wind turbine and characterize the different 
working regions observed in power-speed plane. After this. 
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FIG. 3: Performance curve for one wind turbine in Pinhal Inte- 
rior, Portugal. Circles show all the historical datapoints used in the 
Markov chain modeling i5j and the boxes display the state discretiza- 
tion. 



the same method when appHed separately to both wind speed 
and power production one is able to derive the performance 
curve. Further, the same analysis also provides insight con- 
cerning specific features of the turbine system studied. Sec- 



tion VI concludes this paper 



II. PROPERTIES AND GENERATION OF THE DATA SETS 



FIG. 2: Time series for (a) the power production P of the wind tur- 
bine, (b) the magnitude of the wind speed v. All data series were 
generated with the Markov chain model O described in Sec.|ll] In- 
sets show each time-series in a short time period. All properties are 
normalized to the observed intervals [0, Pmax] and [0, Vmax] respec- 
tively. 



we test our approach to uncover the functional dependence 
of the well-known performance curve, which describes the 
functional dependence of the power production with the wind 
speed. Whereas previous reports have pointed out the benefits 
of deriving the power curve from the drift field||9l, we addi- 
tionally take the diffusion field into account and find that this 
procedure creates additional insight. 

We use the data sets generated by the Markov chain model 
described by Lopes et al. 15]. Using such synthetic data sets 
that properly reproduce the statistical features of empirical 
data sets, allows us to use data sets as large as needed for 
our analysis. Moreover, the Markov chain model behaves as a 
sort of filter for our modeling. 

We start in Sec. |ll]by describing the empirical data used 
to define the Markov chain model as well as the data gener- 
ated with it. In Sec.lnilwe describe our stochastic method for 



analyzing the data and in Sec. IV we apply it to analyze the 
performance curve of a wind turbine. In Sec. [V] we show that 



The data analyzed in this manuscript was simulated from a 
set of measurements from a wind turbine in Portugal, region 
of Pinhal Interior. The measured properties were the power 
production P of the wind turbine and the wind speed v. The 
wind turbine was selected in an eolic park out of a total of 57 
wind turbines. Figure [T] shows an overview of this Portuguese 
eolic park. The time increment between two successive mea- 
sures was At = 10 minutes and the time period covered starts 
in January 1st 2009 and ends in December 31st 2010, yield- 
ing approx. 10^ data points. It has to be remarked that this 
measurements are acquired directly from the top of wind tur- 
bine (nacelle) and might not be optimal for the reconstruction 
of the underlying physical processes for two reasons: first, 
the wind speed measurement is acquired at a point located 
downstream of the turbine blades and can neither account for 
the spatial extension and inhomogeneity of the wind field nor 
for its complex aerodynamical interaction with the turbine 
blades [lOJ. Secondly, the 10 minute sampling period of the 
historical dataset does not allows to resolve the time scales 
of either the turbulent interaction between wind and turbine, 
nor the quick action of the controller system response. Fi- 
nally, missing data records, a low number of data points ifTTI 
or large sampling intervals fl^, and periodicities due to the 
daily and seasonal variations in wind flow often hinders a di- 
rect stochastic analysis of these datasets. Thus, the challenge 
is to devise new methods that can make most use of the in- 
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formation present on these datasets — given that most of the 
data acquisition systems on existent wind farms are Umited — 
with the aim of understanding the dynamic processes of the 
wind power generation, which hopefully can lead to econom- 
ical benefits (production scheduling and maintenance). 

To overcome some of these problems, we employ a re- 
construction of the original data set through a Markov chain 
model. The Markov chain model is established using a joint 
discretization of the wind speed and power variables for the 
state definition and a multiple-step maximum likelihood esti- 
mator for the determination of the transition probabilities |5|. 
In our specific case within the range of each variable P and 
ti we select 80 and 60 states (values). States without a real- 
ization in the time series are deleted. Figure [3] shows the dis- 
cretization of the dataset where the combined states for power 
and speed are indicated with boxes and circles represent the 
historical dataset points projected into the speed-power plane. 

Based on the Markov chain transition matrix P with 
P(i, j) = pij being the probability of transition from state 
to state Sj, the synthetic data sets were generated using the fol- 
lowing Monte Carlo approach. We find the cumulative prob- 
ability transition matrix P^um with Pc„„i(i, j) = J2i^iPt,k 
and select randomly an initial state s^. A random number e 
between zero and one is then uniformly selected and a new 
state Snew is chosen such that Pcumih new) > e. For details 
see Ref. |6|. Figure [2] shows the generated time series. 

The resulting synthetic data series for power production and 
wind speed preserve their cumulative distributions and persis- 
tence statistics, namely the average duration of power produc- 
tion on a certain level (mean sojourn times). Power production 
and wind speed are presented as fractions of the maximum 
observed power Pmax and wind speed Vmax respectively, as- 
suming therefore values between zero and one. 

The synthetic datasets were generated with 2 x 10^ data 
points and show stationary behavior, i.e. have constant moving 
averages (not shown). 

This approach has several advantages over the direct anal- 
ysis of historical datasets. First, high-quality data series of ar- 
bitrary length can be generated, which increases the accuracy 
of the Markov analysis. Secondly, the generated data are by 
construction Markovian, with the reconstruction through the 
Markov chain acting as a filter that removes both noise cor- 
relations and periodicities. Finally, non-Gaussian transition 
probabilities between the states are preserved, which enables 
to study them through higher Kramers-Moyal coefficients. 



III. STOCHASTIC ANALYSIS OF WIND TURBINES 

The co-evolution of two or more stochastic variables, such 
as wind speed and power production can be described through 
a system of coupled stochastic equations, each one defined by 
a deterministic contribution - the so-called drift - and stochas- 
tic fluctuations from possible stochastic sources. In this sec- 
tion we present the general framework to analyze our data and 
in the next section we apply it to the power production and 
wind speed variables. 

For the general case of K stochastic variables, Xi, . . . , Xk 



the vector X{t) = {Xi{t), ....Xxit)) defines the state of the 
system under study at each time instant t. The evolution of the 
state vector yields a stochastic trajectory in phase space and is 
given by the so-called Ito-Langevin eguationllTl [131 [141 : 



~dt 



h(x) + g(x)r(0, 



(1) 



where r(t) = {Ti{t), . . . ,TK{t)) is a set of K indepen- 
dent stochastic forces with Gaussian distribution fulfilling 
the following conditions: {Ti{t)) = and {ri{t)rj{t')) = 
25ij6{t ~ t'). Function h = {/i^} in Eq. (jl]) is the determin- 
istic contribution, describing the physical forces which drive 
the system, while g = {.<7ij} describes the amplitude of the 
stochastic sources of fluctuations r flSl . 

The evolution of the stochastic variables in time yields a 
joint probability density function (PDF), /(X), that evolves 
according to the so-called Fokker-Planck equation 



dt 



i:|-[a'"(x)/(x,o 
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A(f(X)/(X,i) 



(2) 



where the functions D^-^^ and D.-?^ are related with the func- 



tions hi and above, namely 

N 

D^fiX) = ^5,,(X)5,,(X) 



(3a) 
(3b) 



and are usually called drift and diffusion functions respec- 
tively. 

Drift and diffusion functions can be directly derived from 
observed or generated data [15, J_6J, and this fact is the heart of 
our framework. Indeed, the drift and diffusion functions of the 
underlying process are defined through conditional moments, 
namely Itl3i : 



D('^)(X) 



1 M('=)(X, At) 

lim 

At^o At kl 



(4) 



where M^*^' are the first and second conditional moments 
(k = 1,2). These conditional moments can be directly 
derived from the measured data as ||T5l[T6l Mj;^\x, At) = 
{Ydt + At)-Y,{t)\Y{t)^X) and M^f(X,At) = 
((r,(t + At) - Ut)){yj{t + At) - Y,{t))\Y{t) - X> 
where Y{t) = {Yi{t), . . . ,YM{t)) is the A^-dimensional 
vector of measured variables and (•|Y(t) = X) symbolizes 
a conditional averaging over the entire measurement period, 
where only measurements with Y(t) = X are taken into 
account. Important conditions to hold are (i) the underlying 
process is stationary and (ii) the Markovian property is 
fulfifled. 
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FIG. 4: The drift and diffusion coefficients defining tiie co-evolution of P and v: (a) hv, (b) hp, (c) Di,v , (d) D^p = Dp^ and (e) Dpp. The 
PDF for both variables is shown in (f). 



Numerically h and G are determined on a ni x ... x nj^ 
mesh of points in phase space, as a function of the vari- 
ables Xi, using the drift and diffusion functions. Locally, at 
each mesh point, one can always diagonalize matrix G(X) 
and compute their K eigenvalues and K eigenvectors. As 
shown previously 121 Us] [19], this analysis provides informa- 
tion about the stochastic forces acting on the system. Namely, 
the eigenvalues indicate the amplitude of the stochastic force 
and the corresponding eigenvector indicates the direction to- 
ward which such force acts. In a previous work|7| we argued 
that to each eigenvector of the diffusion matrix one can asso- 
ciate one independent source of stochastic forcing Ti and thus 
the eigenvectors can be regarded as defining principal axes for 
stochastic dynamics. In particular, if one eigenvalue is very 
small compared to all the others, the corresponding stochastic 
force can be neglected. In the following sections we present 
a different implication of this principal stochastic component 
analysis, which emphasizes that the vanishing of one stochas- 
tic direction is in fact an indication of a strong functional de- 
pendence between the pair of variables being analyzed. 



following equations: 

dv 



IV. THE PERFORMANCE CURVE FOR WIND TURBINES 



In this section we focus on two variables solely, power pro- 
duction P and wind speed v. Since both series are stationary 
and Markovian we assume them to evolve according to the 



dt 
dP 

lit 



K{v, P) + P)ri + g,p{v, P)T2 (5a) 

hp{v, P) + gp,{v, P)ri + gpp{v, P)r2.(5b) 



In general, the six functions defining vector h and matrix g 
depend on both variables and describe the coupling between 
each other Based on Eqs. p), we derive both h and g from 
the functions D^^^ and D*^^ which, in turn, were extracted 
directly from the synthetic data-set by computing the corre- 
sponding conditional moments using Eq. Q. Note that, solv- 
ing Eq. ( 3b I for computing the matrix g yields multiple solu- 
tions. If g is a solution then all matrices of the form g = gO 
where O is an orthogonal matrix (OO-'^ = 1) are also admis- 
sible solutions. 

In our analysis we take g as the "square root" of matrix 
D(^\ i.e. we diagonalize D^^^ through a proper permutation 
matrix and since all eigenvalues are positive (D*^^) is positive 
definite) we take the square root of each eigenvalue and trans- 
form the matrix back. 

Figures |4^-e show the five components of D*^^^ and D*^^'', 
i.e. the numeric results for both the drift and the diffusion co- 
efficients computed directly from the generated P and v time- 
series. The large fluctuations in the region near to maximum 
power production and wind speed are due to lack of observa- 
tions. Indeed, the joint PDF for P and v (Fig. |4p shows that 
this region is poorly sampled. 

To extract valuable information, next we treat these func- 
tions separately. Namely, we consider the drift vector field 
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FIG. 5: (a) The drift vector D^^' (P, v) = {hp{P, v), h^{P, v)) (see Eq. ijs}) in each (P, v) box used to generate the data (see Fig. [3}. Three 
regions can be identified: Region I having slow dynamics and regions II and III with fast dynamics (see text). Interestingly, the fixed point spot 
in Region II coincides with the wind turbine rated speed (see text). The black curve indicates the performance curve and crosses the bins for 
which D'^' {v, P) vanishes. In the inset, the marginal probability density function of the wind speed v is well fitted by a WeibuU distribution 
with scale parameter A ~ 0.25 and shape parameter k ~ 1.66 and a mode Vmode ~ 0.14wmaa;l iV|. Velocities above Vth ~ 0-52vmax are 
rarely observed, which explain the observed drifts in Region I (see text), (b) Diffusion ellipses in the power production and wind speed state 
space, plotted together with the distribution of data pairs along the boxes shown in Fig. |3] At each box center, the corresponding diffusion 
ellipse is defined by the two orthogonal eigenvector of the diffusion matrix D*^^' computed at that center. The principal axis defining the ellipse 
are aligned along the eigenvectors with a length proportional to the corresponding eigenvalue. Along the performance curve which gives the 
functional dependence between both variables, P and v, the diffusion ellipses degenerate to a line segment tangent to the curve at each box 
center. This feature enables one to use the diffusion matrix of any set of variables for deriving their functional relationships (see text). In the 
inset one sees the ratio of both eigenvalues Xmin/ ^max, using a gray scale (0.2 for black, for white). 



{hy, hp) and the eigenvectors of the diffusion matrix associ- 
ated to its eigenvalues Xmax and ^min- Figure shows the 
drift vector field in the power production and wind speed state 
space, restricted to the sampled region defined by the power 
production curve in Fig. [3] The solid black line is the perfor- 
mance curve computed from {P, v) the joint probability den- 
sity function, shown in Fig. and defines the most likely 
power production for a given wind speed. Three different re- 
gions can be identified. 

Region I is characterized by a large wind speed, i.e. above 
a "threshold" velocity vth that exceeds the turbine rated wind 
speed. Here, the expected behavior of the wind turbine is to 
maintain the power production since there is an surplus of en- 
ergy in the airflow. In this region, the performance curve is 
roughly constant at ^ 0.95Pmax- Still, positive power drifts 
are observed whenever the power production is below the per- 
formance curve. The wind speed drift is large in magnitude 
and always negative, i.e. the drift points towards lower wind 
velocities. 

Region U is characterized by production levels above 
O.bPjnax and wind speeds in OAvmax < w < Vth- Here, there 
is an broad area around the performance curve where the drift 
is relatively small in magnitude. An interesting feature is the 
location of the encircled area in Fig. [5^ near the vth where 
there is a change in the drift vector consistent with an attrac- 
tion point (sink). This shows that the wind turbine rated speed 



was coiTectly selected for this location. 

Region III is characterized by a power production below 

0. BPmax and for wind speed below OAvmax, where there is 
another attraction point at v„iode ~ OAivmax- The separation 
between Region II and III coincides approximately with the 
boundary of the basins of attraction of each stable fixed point, 

1. e. the drift field in Region II and III tends to point away from 
the boundary towards the corresponding stable fixed point 

In previous works |@] the drift vector field around the per- 
formance curve is parallel to the power production axis. In 
Fig. |5^ the vector fields tends to be inclined near the perfor- 
mance curve, because the data analyzed was sampled with a 
much smaller frequency, and therefore the time between suc- 
cessive measures is sufficiently large to observe the conver- 
gence to the stable fixed points. 

Such observations can be more clearly understood by con- 
sidering Fig. together with the marginal PDF of the wind 
speed shown in its inset. The distribution of observed values 
for the wind speed typically follows a Weibull distribution, as 
is known from the literature| 17 1. 

The two dominant trends identified in region I are compat- 
ible with the expected behavior of the power production con- 
trol system present on the wind turbine. For high values of 
wind speed, the controller action upon the blade aerodynam- 
ics is capable of sustaining the production level despite the 



expected decrease of the wind speed. One has to consider, 
however, the time scales involved. With a lOmin resolution of 
the original data, it is not possible to directly observe the rapid 
controller action on the blades, only the controller actions that 
occur on larger time scales, such as the rotation of the tower. 
However, even at large samphng times, the data set catches 
some events outside of the power curve and the subsequent 
conditional moments mirror the controller action that forces 
the system back on the curve. 

A closer look in Region II enables one to identify a fixed 
point region (D'^^ ~ 0) at high power production levels and 
wind speed v ~ O.SVmax- This speed value coincides with 
the rated wind speed, i.e. the speed for which the turbine was 
designed and at which it operates at an optimal regime. It 
can therefore be concluded from our analysis that the turbine 
has been well adjusted, and it remains to be seen if similar 
conclusions can be drawn when applying our method to arrays 
of turbines. 

Another important application of our method deals with the 
diffusion matrix. As explained in the previous section, by di- 
agonalizing the diffusion matrix at each point of the phase 
space one is able to determine the two eigendirections for dif- 
fusion. Being orthogonal to each other, these two eigendirec- 
tions define an ellipse with major and minor axis proportional 
to the corresponding eigenvalue. Figure [5J5 shows the diffu- 
sion ellipses in phase space. Region I is characterized by the 
largest ellipses indicating very large fluctuations, while there 
is an area in region II, which present small fluctuations and 
corresponds to the fixed point areas identified in the drift. In 
the high-slope region of the power curve, the ellipses degen- 
erate, i.e. one eigenvalue is negligible when compared to the 
other {Xmin/^max ^ 0). The inset of Fig.jsJ? shows in a gray 
scale the quotient Xmin/Knax between the smallest and the 
largest eigenvalue. White corresponds to zero quotient, while 
values in [0.2 . . . 1] are colored in black. Clearly, a white curve 
can be identified in the inset which follows the performance 
curve shown in Fig. [3] 



V. DERIVING THE PERFORMANCE CURVE FROM 
UNIVARIATE STOCHASTIC DYNAMICS 

The inset of Fig.jSj) shows that along the performance curve 
one eigenvalue is typically much larger than the other As we 
will argue next, such feature is an indication that in fact P is a 
function of v, which in the case of power production and wind 
speed yields the performance curve drawn in Fig. [3] and [5] 

Interestingly, we will see that analyzing both wind speed 
and power production separately enables one to extract valu- 
able insight about the full dynamics and behavior of the wind 
turbine with the atmospheric wind. In fact, the 2D analysis of 
the performance curve summarized in Fig.[5]can indeed be ac- 
cessed through a one dimensional stochastic analysis of each 
variable P and v separately. 

To see this one first takes P and v as two general variables 
fulfilling Eqs. (jsjl and observes that if P = P{v) there are not 
two independent stochastic forces, but only one, yielding for 
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v/v , P/P 

max max 

FIG. 6: Uncovering properties of wind turbines by analyzing data 
series of wind speed v and power P separately (see text and Eqs l|6j 
and lj7|). Horizontal axis indicates the value of v/vmax for hi, and 
and the value of P/Pmax for hp and gp. 

dv ~ 

— = hy{v) +gy{v)T, (6) 

and for P 

dP 

— ^hp{P) + ~gp{P)T, (7) 

where functions h and g are of course different from the drift 
and diffusion functions defined above in Eqs. (jSjl, since only 
one variable is taken into consideration for the stochastic mo- 
tion equation. 

Figure |6] shows the drift and diffusion of both the wind 
speed and of power production determined for the model de- 
scribed by equations Eqs. The drift of the wind speed, 
has three zeros. These zeros correspond to three fixed 
points, two stable {v ~ v^node ^ 0.14w„ax and v ~ vth), 
and one unstable at u ~ QAvmax- Thus, for wind speed 
below ~ QAvmax, the airflow is unstable and unsuited for 
power production, while wind speeds above ~ QAvmax pro- 
mote power production. The first zero of hy indicates ap- 
proximately the mode of the wind speed value distribution - 
compare with inset in Fig. [5^ - and the other two zeros mark 
the transition between two different regions identified above 
in Fig. |5] The transition between Region I and Region II is 
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FIG. 7: Obtention of power production P as function of wind speed 
V by integration (see Eq. l[8|). Both Eqs. j|9]( are fulfilled, analyzing 
both series P and v separately (see Eqs. ^ and jTj). 



marked by vth- The transition between Region 11 and III is 



more subtle and deals with the zero at OAvj. 



and with the 



functions for both variables: it is located at the value Pi, which 
corresponds to v ~ QAvmax, cf. Fig. [5^. Moreover, Fig. [6^ 
also shows that positive drifts are located for small wind speed 
(up to Vmode) and for Region II (see Fig. [5]l, while Fig. |6|3 
shows the transition value Pf, of power production is close to 
the maximum of gp. 

The transition between Region III and Region II is located 
at a minimum of the diffusion for the wind velocity. As 
shown in Fig. [S]?, here the drift changes to a positive value. 
In other words, above v ~ QAv,nax the expected change of 
the wind speed is towards higher values. The drift reaches a 
maximum on Region II prior to a steep change towards neg- 
ative values. In this region, the wind speed values are not in 



the range of the extreme weather conditions and are also not 
as frequent as the lower wind speed values. However, it is 
frequent enough to be associated with a commonly repeated 
pattern, i.e. the daily pattern seasonality, in which the air- 
flow is mainly induced by thermal differences. This pattern is 
responsible for most of the power production in this wind tur- 
bine and the main reason for a second attraction point at Vth- 
For higher wind speed values, the drift changes again to neg- 
ative values (Region I) since very high wind speed is usually 
of short duration, i.e. extreme wind gusts. 

Bringing all the above observations into account, one con- 
cludes that there is a strong agreement between the regions 
defined in the context of Fig. |5] and the sign of the wind speed 
drift. 

Having analyzed separately both properties, v and P, we 
end this section by showing that from the drift and diffusion 
coefficients, h^, g^, hp and gp, one indeed obtains the func- 
tional dependence between power production and wind speed. 
To that end we assume that Eq. (j6]l holds for v and that the 
other variable P is an exclusive function P{v) of v. Thus, we 
can take the Ito-Taylor expansioni 14 1 of its differential 



dP{v) 



P{v + dv) - P{v) 

—dv + ^-r^dv 
av dv^ 
dP~ 



d^P 



dv 



dv'^ 



9v 



0{dv^ 
dP 



dt 



dv 



g^dw (8) 



using the differential dv — hy{v)dt + gy{v)dw. Therefore, 
identifying 



tip — ——tiy 
dv 



d'^P 



dv 



2 



9P 



dP, 
dv 



gv 



(9a) 
(9b) 



which can be solved with respect to the two derivatives of P{v) yielding the numerical integration scheme as follows: 

P{v + At) 



dP ^ d'^p 

P{v) + ^|.,P(.)A« + -^|„,p(„)(At;)2 + 0{{/\vf) 

1 hp{P{v))~gy{v) - K{v)~gp{P{v)) 



9v{v) 



(Aw)2 ^-©((Ai;) 



(10) 



Figure [t] shows the integration of dP{v) for the condition 
-fo(''^o) = for vq = 0. The deviations can be attributed 
to the fact that Eqs. (j6| and (j7]i are strictly only valid in those 
regions where the eigenvalues of the diffusion matrix show a 
large difference between them, Xmin ^ ^max^, cf. inset of 
Fig |5] where 0.2 < i; < vth- In addition, it also holds only on 
the performance curve, and applying Eqs. (j6]) and (j7]i there- 
fore also neglects the asymmetry of the drift functions with 
respect to this curve. 



VI. DISCUSSION AND CONCLUSIONS 



Investigating a wind turbine from a commercial wind park, 
we report the reconstruction of the stochastic performance 
curve in the variables wind speed and power production, us- 
ing both drift and diffusion coefficients. These coefficients, 
describing the respective deterministic and stochastic interac- 
tions of wind field, turbine aerodynamics, and controller ac- 
tion, are estimated from a synthetic time series reconstructed 
from a Markov Chain model of the original measurement data. 
We argue that this reconstruction is highly superior to a direct 
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evaluation of the measurements. 

As a main finding, we present the fact that the reconstruc- 
tion of the power curve using both drift and diffusion coeffi- 
cients uncovers additional information not visible in an anal- 
ysis of the drift field alone|9|, even though we are using mea- 
sured data of a very low measurement rate as model input. 
Specifically, our analysis reveals the existence of various dis- 
tinct regions in the wind speed-power production plane 

In addition, we have been able to reconstruct the power 
curve from the drift and diffusion coefficients, using a method 
which should be able to uncover functional relationships be- 
tween stochastic variables in a wide range of experimental se- 
tups. 

We our convinced that our approach is helpful in obtaining 
a better understanding of the complex dynamics that deter- 



mines power production in wind turbines, and that it can be 
useful in optimizing the output. It is planned to extend this 
methodology in order to consider coupled systems of wind 
turbines in nearby locations. 
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